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Abstract 



& ', We develop a characterization method of electroconvection structures in a 

planar nematic liquid crystal layer by a study of the electric current transport. 



Because the applied potential difference has a sinusoidal time dependence, we 
define two electric Nusselt numbers corresponding to the in-phase and out- 
of-phase components of the current. These Nusselt numbers are predicted 
theoretically using a weakly nonlinear analysis of the standard model. Our 
measurements of the electric current confirm that both numbers vary linearly 
with the distance from onset until the occurence of secondary instabilities; 
these instabilities also have a distinct Nusselt number signature. A system- 
atic comparison between our theoretical and experimental results, using no 
adjusted parameters, demonstrates reasonable agreement. This represents a 
quantitative test of the standard model completely independent from tradi- 
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tional, optical techniques of studying electroconvection. 
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Although spontaneous pattern formation within structureless environments pervades na- 
ture a comprehensive understanding of this complex behavior remains elusive. Therefore 
well-controlled experimental systems exhibiting pattern formation are extensively studied; 
among these, thermoconvection of a layer of fluid heated from below || and electrocon- 
vection of a nematic liquid crystal layer are particularly interesting since they allow 
very large aspect ratio geometries. In both of these systems, convection structures form 
spontaneously when the applied stress, i.e. the gradient of either the temperature or the 
electric potential, exceeds a critical value. These inherently non-equilibrium structures can 
only persist when there is an energy source to overcome the dissipation associated with the 
flow. Therefore, energy transport studies represent a particularly valuable technique for elu- 
cidating the essential nature of the instabilities that lead to these patterns. For example, the 
first accurate determination of the stress necessary to induce thermoconvection was made 
by measuring the heating power required to sustain a desired temperature difference across 
a thin layer of water ||. This power is customarily expressed as the Nusselt number, defined 
as the heat flow across a fluid layer relative to the heat flow required in the absence of fluid 
flow. Nusselt number measurements remain a method of choice for studies not only of the 
structured states that occur during thermoconvection when the stress is only slightly above 
its critical value || but also of the turbulent flow that occurs when the stress is enormous 
0. By contrast the electroconvection of a planar nematic liquid crystal layer, which rep- 
resents a similar but fully anisotropic model pattern forming system, has previously been 
studied only with qualitative or semi-quantitative optical techniques. Reports of energy flow 
measurement during electroconvection are rare ||, and no theoretical studies of the energy 
transport exist for this system. The aim of this Letter is to fill this gap. 

Electroconvection is obtained when an a.c. electric potential, y/2V cos(u?t), is applied to 
two horizontal (_L z) electrodes separated by distance d confining a nematic liquid crystal. 
Here we focus on the planar anchoring case where the director field n is fixed to x at the 
confining electrode plates. The instability relies on a coupling between n, the velocity field, 
v and the induced charge density, p e or equivalently the induced electric potential <f), such 
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that the full electric field reads E = \plV jd [cos(ut)z — dV0]. At moderate frequencies u, 
when V exceeds a critical value V c , the instability sets in the form of normal conduction 
rolls of wavevector q = gx ; at large frequencies dielectric rolls are observed but we do 
not consider this regime in this work. These phenomena are well explained via the standard 
model (SM) for electroconvection HQ where the charge conduction in the liquid crystal is 
assumed Ohmic. In addition to linear properties (values of V c and q as a function of uj), 
the SM explains several secondary instabilities that are experimentally observed, such as 
the transitions to zig-zag rolls, stationary and oscillatory bimodal patterns and abnormal 
rolls || . The one phenomenon which the SM has been unable to predict is the traveling roll 
state. For this, a new approach, the weak electrolyte model, was developed, in which the 
electrical conductivity is assumed to be due to two species of dissociated ions having different 
mobilities [TIJ. With this model a semi-quantitative agreement with experimental results 
on traveling rolls was demonstrated |T0| but fitting parameters were necessary. Because 
traveling rolls are not encountered at low frequencies where we performed our experiments, 
and because the SM is much less complicated, we did not use the weak electrolyte model. 

The total current / through the nematic cell enclosed by the horizontal electrodes of 
area S can be calculated as the circulation of the magnetic induction H. From the Maxwell- 
Ampere equation, V x H = j + <9 4 D, I is sum of the conduction and displacment currents: 



where, within the SM, j = crj_E + (T a (n-E)n-|-p e v ; D = e_i_E + e a (n-E)n ; a± (cry) are the 
conductivities perpendicular (parallel) to n , and, e± (ey) are the dielectric permittivities 
perpendicular (parallel) to n , a a = cry — crj_ and e a = ey — e± are the corresponding 



{—d/2 < z < d/2) chosen, because of the Maxwell- Ampere equation. In the quiescent (no 
convection) state n = x, v = and = 0, therefore 




(1) 



anisotropies. Note that the surface integral in eq. ([I]) does not depend on the z-value 



1 = 1° = J r ° cos(ut) - J° sin(urt) 



V2VS 
d 



[cr_|_ cos(ujt) — e±u> sin(u;t)] . 



(2) 
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In the convecting state all fields n, v and are modified, as is /. Within the SM, for 
homogeneous stationary roll solutions 

I = I r cos(uit) — Ii sin(u;t) + higher temporal harmonics (3) 

where the amplitudes of the higher temporal harmonics are expected to be much smaller 
than I r and Ii , at least at intermediate frequencies uo ~ 1/tq where tq = e±/a± is the 
charge- diffusion time. We define the real and imaginary reduced Nusselt numbers, M r and 
Mi , as I r /Ir — 1 an d h/lf — 1 , respectively. Thus M r = Mi = in the quiescent state, while 
in the convecting state M r measures the excess energy dissipation due to convection of the 
nematic liquid crystal, that is the time average 

In heuristic terms the effective resistance of the nematic layer is changed by convection from 
Ro = d/(a±S) to R = Rq/(1 + M r ) ; equivalently the imaginary Nusselt number measures 
the change in the effective capacitance of the nematic layer, Co = e±S/d in the quiescent 
state, C = C (l + Mi) in the convecting state. When the reduced distance from onset 
e = V 2 /V 2 — 1 is small, the electric Nusselt numbers can be calculated for homogeneous 
rolls using weakly nonlinear methods. Assuming that the leading convection amplitude, 
A, associated with the linear roll mode, remains small, a systematic expansion in powers 
of A is performed. After adiabatic elimination of the slave modes and calculation of the 
resonant saturating cubic terms approximate roll solutions are obtained together with the 
relation A(e) = aoy^ • The current can then be calculated from eq. ([I]). For symmetry 
reasons the first contribution from the convection modes comes at order A 2 , and therefore 
one expects M r ~ A 2 ~ e , Mi ~ A 2 ~ e in the weakly nonlinear regime. That is, 
Nusselt numbers allow a direct measurement of the convection amplitude A, and therefore 
a test of the supercritical law A(e) = a^^fe . In order to make this clear and to obtain 
approximate analytic formulae, one can first use the quasi-unidimensional approximation, 
where all fields are considered at the middle of the layer (z = 0) and only their x-dependence 
is kept. The linear normal roll mode then assumes the form n z = —AN z sin(qx),v z = 
A/(qr )V z cos(gx), (p = Aj \qd)[<& c cos (out) + $ s sin(a;t)] cos(ga;) where, as in the rest of our 
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theoretical calculations, we only keep the lowest nontrivial time-mode for each field. We 
also choose as a normalization condition N z = 1 ; then V z , $ c and $ s are calculated at fixed 
frequency by solving the linear neutral eigenproblem. With p e = V • D, j z + dtD z can be 
easily calculated at z — 0. Keeping only the horizontally homogeneous terms because of the 
surface integral in eq. ([I]), one obtains to lowest order in A 

K = W a N z (N z - $ c ) + e\^ c V z - e' a N z (V z + lut $ s )} (4) 



2 
A 2 



: || 

e' a N z (N z - $ c ) + ^ ( a ' a N z - e\V z [ 

UJT " 
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with o' a = a a /ax_ , e' a = e a /e± and eji = 6||/ej_ . For standard nematic materials with large 
positive a' a (see e.g. eq. (||)), the leading term in M r eq. ([|) is the anisotropic conduction 
term in u' a N^, which imposes a positive value of Af r . One thus expects that the tilt of 
the director out of the plane in roll structures will enhance the electrical conduction of the 
layer and finally the in-phase current. Concerning Mi eq. (|5|) it should be noted that Q s /uj 
tends to a finite positive value when u — > 0. Two terms of eq. (|5|) control the sign of 
Mi . The first term in e' a N z reveals a diminution of the effective capacitance of the cell 
due to the director tilt, since the dielectric anisotropy e' a of the nematic materials used in 
electro convection is usually negative. The other important contribution is the positive term 
in a' a N z Q s /uj , which expresses that the potential modulation induced by the convection 
creates by coupling with the director tilt an out-of-phase current U > (see eq. @). Since 
$ s /u; decreases with u, this positive term can compensate the negative term in e' a N z at low 
frequencies only, and for nematics with large a' a . A numerical calculation, using a standard 
Galerkin technique to expand the z-dependence of all fields in test functions, can provide 
a more accurate evaluation of M r and Mi ■ For this purpose we have modified the code 
developed in |J to use Tchebyshev polynomials as the test functions in order to accelerate 
the convergence (typically 4 z-modes were sufficient), and inserted a procedure to calculate 
the Nusselt numbers. For convenience the current eq. ([!]) is evaluated at the lower plate 
z = —d/2 where, because of the boundary conditions, j z + d t D z reduces to (a± + e±d t )E z . 
Thus, since the convection induced potential is even under z t— > —z at linear order, but odd 



at quadratic order, one sees that the leading contribution to / comes from the potential 
part of the homogeneous quadratic slave mode noted A 2 V 2 (q, — q) in eq. (27) of ||. Of 
course the saturation at cubic order needs also to be calculated in order to provide the law 
A = a$\[i. We will return to the numerical results (Fig. |^), which confirm the trends found 
from the analytic formulae eqs. (Q), (Q), after presenting our experimental results. 

We use the "classical" arrangement based on a pre-fabricated liquid crystal cell Jl2| . 



The glass plates have no spacers, glue, etc. within the active area; they are separated by 
d = 23.4 ± 0.5/xm. With air only between the plates, we measure, using a auto-balancing 
1kHz bridge, the capacitance of the cell in order to determine accurately (within 8 ppm) the 
ratio S/d (nominally 5* = 5mm x 5mm). After this measurement the nematic liquid crystal 
methoxy-benzylidene butyl-aniline (MBBA) ||13|| , doped with 0.0005% tetrabutyl ammonium 
bromide, is introduced between the transparent conducting electrodes. The filled cell is 
placed in a temperature controlled housing, and then introduced between the pole faces of 
a large electromagnet. As the nematic liquid crystal undergoes the magnetically induced 
splay Frederiks transition, the capacitance and conductance of the cell are monitored. ^From 
these measurements we obtain both electric conductivities and both dielectric constants |L4| . 



Hence, we measure in situ all the electrical transport properties of the specific nematic liquid 
crystal used. For the experiments reported here, all at 28°C, we find 

cr ± = (8.5 ± 0.8) lO^lWT 1 , a' a = a a /a ± = 0.35 ± 0.04 , 

e_L = (4.65 ± 0.03)e , e' a = e a /e± = -0.080 ± 0.001 , (6) 

with eo the vacuum dielectric permittivity. After these measurements, the nematic cell is 
transferred to the stage of a polarizing microscope so that shadowgraph |1| images can 
be obtained concomitantly with the electric current measurements. A function generator is 
used to produce a sinusoidal voltage signal which is in turn amplified, and applied to the 
cell. The path to ground for the current traversing the cell is through a current-to-voltage 
converter. The output signal from this converter is measured by a lock-in amplifier, whose 
reference signal is supplied by the original function generator. Before any measurements are 



taken, the nematic cell is replaced by a purely resistive load and the phase setting on the 
lock-in is adjusted to zero the out-of-phase current component. The nematic cell is then 
re-inserted. Then, at a selected frequency V is raised in small steps. At each step, after 



waiting several seconds, I r and U are recorded |16| . This proceeds until a maximum desired 
V (well above the threshold value V c ) is reached. Then, the process is reversed, and the 
currents recorded as V decreased. The difference in current for increasing vs decreasing V is 
less than 2%. When V is raised above V c , the electric current traversing the liquid crystal 
measurably deviates from its value in the quiescent state, 1°. In order to determine V c from 
either the in-phase or out-of-phase current data, we first determine a baseline for If (If) by 
fitting a straight line to I r (ij) vs V for V much smaller than V c ; see Fig. [1} These values 
of If/V and If /V provide independent measurement of cr_i_ and e± (see eq. (||)) that agree 
within 5% with the direct measurement of these parameters using the Frederiks transition. 
The Nusselt numbers as functions of V are then calculated by subtracting unity from the 
ratios I r jlf and h/If ■ By fitting another straight line to M r (Mi) in the region where it 
deviates from zero, we define V c as where this line crosses zero (see the insets in Figs. [TJ and 
0). For one ramp of V, the three values of V c determined from the Nusselt numbers and the 
traditional shadowgraph technique agree with each other within 0.01%. 

Our apparatus did not reach sufficiently large V to measure the crossover to the dielectric 
regime (see e.g. f§]); we therefore estimated a characteristic "cutoff" frequency uj c by fitting 
V c (uj) to the function A/(uj — u> c ). We found typically u c /(2ti) = 645Hz, but during the 
course of taking the measurements (2-3 months) this quantity varied by ±10%, probably in 
connection with variations of the electrical parameters, especially the conductivities. 

In Fig. we plot both the real and imaginary Nusselt numbers vs e. While M r is always 
observed to be positive, for the data set shown Mi is negative. In some cases (discussed 
subsequently) it becomes positive. Note also that M r is at least ten times larger in magnitude 
than Mi- 
Close to threshold, i.e. for < e < 0.1, both Nusselt numbers are proportional to 
e as shown in the inset for the real Nusselt number: this confirms the supercritical law 
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A ~ y/e. The variations of the corresponding slopes M r /e and Mi/e vs ut are given on 
Fig. which represents the results of several ramps in e at each frequency. For e ^ 0.1 
the curves deviate from straight lines; the "knees" in the curves indicate clearly the onset 
of secondary instabilities. Specifically, the two arrows shown on Fig. |2| correspond to the 
onset potential differences for the zig-zag instability []I7[ and the spontaneous generation of 



dislocations associated with the so-called "defect chaos" |18[ . 

To compare our experimental results with the SM calculations, we used the elastic con- 
stants and the viscosities measured for MBBA at 28°C in and the electric parameters 
that we determined independently (eq. (§)). Varying those parameters within the stated un- 
certainties, we calculate (with the fitting procedure defined above) for the "cutoff" frequency 
u c / (2tt) = 730 ± 120Hz, in agreement with the measured value. Systematically varying the 
electrical parameters within the experimental error bars, we also calculate, with the weakly 
nonlinear numerical code introduced above, the bands of possible values of M r /e and Mi/e. 
These bands are drawn in gray on Fig. []. The extremal values turn out to be obtained 
by variation of only a' a , with the upper (lower) curves for both M r /e and Mi/e obtained for 
the largest (smallest) value of o' a . This is consistent with the fact that the leading positive 
terms controlling the Nusselt numbers as seen in the approximate formulae eqs. (^), ([S]) 
are proportional to a' a . Note also that for large a' a we expect M% to be positive at low 
frequencies, while Mi is always negative for a' a small . There is a good agreement between 
experiments and theory concerning the imaginary Nusselt number Mi ; on the other hand 
the real Nusselt number M r decreases more abruptly in the experiments than in the theory. 
However the agreement obtained for Mi at all frequencies and M r at small frequencies is par- 
ticularly significant since (contrarily to the standard approach in nematic electroconvection 
where usually a± is fitted) no adjusted parameters have been used. 

In conclusion, electric Nusselt number measurements are validated as a new and powerful 
method of characterization of electroconvection. This wholly quantitative technique stands 
in contrast to traditional optical methods which only become quantitative in certain limit 
cases. This technique affords a precise determination of the threshold voltage for the onset 



of electroconvection as well as secondary instabilities. Nusselt number measurements also 
represent an important quantitative tool for testing competing theoretical descriptions of 
electroconvection. Here, with only the limitation of relying on tabulated values of some 
material parameters, we have shown that the standard model for electroconvection gives 
satisfactory predictions of the Nusselt numbers near onset. One conspicuous explanation for 
the remaining discrepancies may be that nematic liquid crystals are quite clearly electrolytic 
conductors, and thus the Ohmic conduction assumed in the standard model introduces an 
important approximation. Thus, it clearly is of interest to extend the calculations presented 



here within the so-called weak-electrolyte model |]10| . Future directions of this work include 
also employing liquid crystal materials for which the applicability of the weak electrolyte 
model has been established, and experiments in the highly nonlinear, dynamical scattering 



regimes that occur at very large e [21 
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FIGURES 

FIG. 1. In-phase current I r vs the applied voltage V at a frequency lo/(2tt) = 100Hz i.e. 
ujtq = 0.30. Inset: blowup of the neighbourhood of V c . 

FIG. 2. Real and imaginary reduced Nusselt numbers vs the distance to threshold e for 
electroconvection of MBBA at the same frequency than Fig. The arrow on the left indicates the 
onset of the secondary zig-zag instability and the arrow on the right indicates the onset of defect 
chaos. Inset: blowup of the e ~ region indicating how M r initially increases linearly with e. Small 
pre-transitional effects caused by sample imperfections are also seen. 

FIG. 3. Black circles: measured ratios M T j e and A/i/e at small e vs the dimensionless frequency. 
The vertical error bars correspond to the imprecision in measuring M je and to the dispersion 
between different e ramps, while the horizontal error bars originate from the variation in the 
charge-diffusion time tq . Gray bands: intervals of values of the same quantities deduced from the 
weakly nonlinear SM calculations, taking into account the variations of the electrical parameters 
of the nematic liquid crystal eq. @. The upper (lower) theoretical curves on both graphs were 
calculated with the largest (smallest) value of a' a = a a /cr±. 
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